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ABSTRACT 

Although coherent large-scale structures such as filaments and walls are apparent 
to the eye in galaxy redshift surveys, they have so far proven difficult to characterize 
with computer algorithms. This paper presents a procedure that uses the eigenval- 
ues and eigenvectors of the Hessian matrix of the galaxy density field to character- 
ize the morphology of large-scale structure. By analysing the smoothed density field 
and its Hessian matrix, we can determine the types of structure - walls, filaments, or 
clumps - that dominate the large-scale distribution of galaxies as a function of scale. 
We have run the algorithm on mock galaxy distributions in a ACDM cosmological 
N-body simulation and the observed galaxy distributions in the Sloan Digital Sky Sur- 
vey. The morphology of structure is similar between the two catalogues, both being 
filament-dominated on 10-20 h Mpc smoothing scales and clump-dominated on 
5 hT x Mpc scales. There is evidence for walls in both distributions, but walls are not 
the dominant structures on scales smaller than ~ 25 hr l Mpc. Analysis of the simula- 
tion suggests that, on a given comoving smoothing scale, structures evolve with time 
from walls to filaments to clumps, where those found on smaller smoothing scales are 
further in this progression at a given time. 



1. Introduction 



Inflationary models of the early universe (iGuthll 19811) suggest that the large-scale distribution 



of matter should be well described by a gaussian random field, the statistical properties of which 
can be completely specified by the power spectrum. As structure evolves, however, non-linear 
effects become increasingly important and higher-order statistics must be computed if we wish 
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our description to be complete. A number o f studies have been performed to measure the three- 
point correlation function of galaxies (e.g. IPeebles & Grothl 1 19751 : iGaztanaga & Friemanl 1 1994 : 

Nichol et al]|2007l : Kulkarni et al.ll2007l) and its Fourier Transform (the bispectrum), as well as N- 

point statistics (e.g jFry & Gaztanagalll993 : Verde. Heavens. & MatarreselboOO : Ross. Brunner. & Myers 
20071) . However, the filaments, walls, and clusters that appear in the non-linear regime produce 
very high-order correlations which are difficult to compute and interpret. 



Three-dimensional maps of the galaxy distribution show walls and filaments (e.g., [T hompson & Gregory 



19961 : lAragon-Calvo et al.ll2007l : 



Colless et al 



GottetaL 



19781 : Ide Lapprent. Geller. & Huchrall 19861 : teller & Huchrall 19891: 

2005b as do cosmological N-Body simulations (e.g ., Davis et al.ll985 : Bond. Kofman. & Pogosyan 



2001; 



Hah n et al.ll2007l). At present, there are a number of methods for 



quant i fying filamentari t y (e.g.. Klypin & Shandarin 
20041 : Ishen et all EoOfjl; iNovikov. Colombi. & Dorel 



1993; IShandarin & Yessl 1 19981 : IStarck et al 



2006), a s well as for tracing individual fil 



aments (e.g. , Sousb ie et al. 120081 : lAragon-Calvo et al.ll2008l : ISousbie. Colombi. & Pichonl 120091 : 



Bond. Strauss. & Cenll2010h . Filaments produce deep potential wells on large scales, surpassed 



only b y clumps^J. As such, they will present a si gnal for gravitational lensing studies on the largest 
scales (IDietrich et al.ll2005l : iMassey et al.ll2007l) . Furthermore, filaments can act as a pathway for 



matte r accreting onto young galaxies and galaxy c lusters (e.g. iTanaka et al.l 120071 : iRoberts et al. 
20071) and can align the spins of dark matter haloes (|Hahn et al.ll2007|) . 



In the classic picture put forward by IZerdovichl (I1970I) . an ellipsoidal overdensity will col- 
lapse along its shortest axis first, forming a 'pancake' -like structure. Some early models of struc- 
ture formation, such as the 'hot dark matter' models, led to a 'top-down' scenario, in which an 
excess of power at large scales led to a universe filled with giant sheets, extending for tens t o 
hundreds of megaparsecs (IZel'dovich. Einasto. & Shandarinlll982l : iFrenk. White. & Davisll 19831) . 
We now think that dark matter is cold and that structure formed in a bottom-up fashion, with the 
smaller structures forming before the larger ones. Walls may still exist, but they are not as obvious 
in either the galaxy catalogues or cosmological simulations as they would be in the hot dark matter 
models. Furthermore, it is not clear that we can apply the simple picture of ellipsoidal collapse to 
structures that are experiencing tidal forces from nearby structures and continual fragmentation on 
smaller scales. 

One approach to desc ribing bottom-up structural evolution is known as the 'peak-patch' the- 
ory (IBond & Myerslll996|) . in which clumps, the rare overdensities in a smoothed density field, 
form first. Filament-like features then arise along the longest principal axis of the tidal tensor, 
forming bridges between the clumps, followed by wall-like structures connecting the bridges. Al- 



Clumps are defined as nearly spherical overdensities that have undergone significant collapse along all three 
principal axes. Bound clusters are a subset of this class of structure. 



-3- 



though this description incorporates tidal interactions between neighboring structures, it focuses on 
only one smoothing scale at a time and is forced to resort to increasingly complex Appoint statistics 
to describe the development of structures on scales much larger than the smoothing length. 

In this paper, we are focused on measuring the dimensionality of structures in the overall 
distribution of matter as a function of scale. Past studies of the correlation function have shown 



large-scale structure to b e consistent with a fractal on < 10 h 1 Mpc scales ([Martinez & Coles 



1994l : IPan&Colesll200Ql) . suggesting that many of the structures apparent to the eye are nested 



into larger structures. Any measure of structure that ignores this fact is giving an incomplete 
picture of the matter distribution. 

In what follows, we will describe a procedure to quantify the prominence of structures, such 
as filaments or walls, in the large-scale distribution of galaxies and then apply this procedure 
both to simulations and to o bservations of th e large-scale distribution of galaxies from the Sloan 
Digital Sky Survey (SPSS: lYorket all boook First, in § [2l we will define the local structural 



parameters. Walls, filaments, and clumps will each have a unique fingerprint in the resulting 
parameter space (A -space) and we will present these fingerprints, along with their dependence on 
the properties of a given structure, in § [3j In § HI we describe the N-body simulations and mock 
galaxy catalogues with which we will develop our statistics and compare the large-scale structure 
data. In Sections [5] and [6l we discuss the A -space distributions of gaussian random fields and dark 
matter in a cosmological simulation, respectively. We will examine these statistics and compare 
these distributions to those of real galaxy data in § [7j Finally, in §[8l we will summarize the results 
and discuss their implications for our understanding of large-scale structure. Bond, Strauss, & Cen 
(Paper II, in preparation) will describe a method of finding individual filamentary structures. 



2. Method 

2.1. Single-scale smoothing of the density field and its Hessian 

Filaments, clusters, and walls all present sharp features in the density field along at least one 
of their principal axes. If we wish to bring out such features, we should generate not only the 
density field, but also its second derivatives. The smoothed density field is defined as, 

p(x) = Jf(x-x') P (x')d 3 x', (1) 

where p(x) is the density field and f(x) is the smoothing kernel. The density field is initially 
composed of a sum of delta functions (at the positions of the galaxies or dark matter particles), 
which are zero-valued at infinity, so integration by parts yields the smoothed Hessian (i.e. matrix 
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of second partial derivatives), 



j f{x-x' 



dx[dx\ 



d 3 x' 



d 2 f(x-x') 
dx[dx\ 



p{x')d 3 x'. 



(2) 



A gaussian smoothing kernel will be used throughout this paper. 



The convolution is performed on a discrete grid after first performing a Fast Fourier Transform 
(order MogAO on the kernel and data vectors. This operation implicitly assumes periodic boundary 
conditions, a criterion that is met by the simulation data (which is in a periodic box), but not by 
the survey data. For the latter case, we will only search for structures more than two smoothing 
lengths from the box edge. 

The idea to use the Hessian matrix to characterize local large-scale structure was pioneered by 
Colombi. Pogosyan. & Souradeed (12001) and later used for the identification of individual large- 
scale structures by Aragon-Calvo et alT ( 2007). A similar method, which uses the eigenvalues of 
the tidal shear field, was implemented by lHahn et al. ( 2007 ) to characterize large-scale structure in 
the vicinity of dark matter halq s in N- body simulations. This technique was later applied to galaxy 
redshift surveys by iLee & Led (|2008|) . The method presented here is also similar to that used by 



Forero-Romero et al 



(120091) to present a dynamical classification of large-scale structures. 



It has been suggested by some authors (e.g. Istein 1997 ) that, because filaments present them- 
selves on a multitude of scales, one should adaptively smooth the galaxy density field. It's true that 
a visual inspection of a high-resolution density map of large-scale structure will reveal filaments 
on scales as small as ~ 1 hr l Mpc, but these filaments will often be embedded in larger structures, 
which will be difficult to identify in a single density map that deconstructs them into their smaller 
components (see § 18. 1 1 for further discussion). For these reasons, we will smooth separately on a 
series of length scales when searching for filaments and walls. 



2.2. The A parameters and axis of structure 

The Hessian matrix can be visualized as an ellipsoid with the major axis aligned along the 
direction of lowest concavity. In filamentary regions, this direction will be along the filament 
itself, so the local orientation of a filament is given by the major axis of the Hessian ellipsoid 
at that point. By diagonalizing the Hessian matrix, we move into the coordinate system of this 
ellipsoid and obtain its principal axes (the Hessian eigenvalues) and orientation (eigenvectors). 

The Hessian eigenvalues at a given grid cell will be denoted Aj and ordered such that X\ < 
A2 < A3. Along with the density and gradient vector, these provide an indication of the shape and 
contrast of local structure. The orientation of the structure at a given grid cell is given by A3, the 
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eigenvector corresponding to A3. Although real filaments will often have a non-negligible gradient 
along their length, we will consider only the axis along the filament, hereafter referred to as the 
'axis of structure' . We wish to keep our notion of a filament very simple - locally, it is a structure 
that is concave down along two principal axes and nearly flat along the other one. We will make 
no assumptions about the sign or magnitude of the first derivative along its length, but for the sake 
of convention, we will define A3 such that 

A 3 -Vp > 0. (3) 
We define the dimensionless eigenvalues, 

K - ^ (4) 

where / is the smoothing length and p is the mean density. Near the centres of filaments, where the 
density field is concave downalong two axes of the filament, we expect the following relationships 
between the A '-parameters: 

X[ <0 

A^ ~ X[ < (5) 

For walls, the density field is concave down along only one axis and slowly varying along the other 
two, 

X[ <0 

I^KIAd (6) 

Filaments are essentially one-dimensional objects, so a single vector can be used to represent 
their orientation. However, if we wish to consider walls, an additional eigenvector, A2, is needed. 
The full set of parameters to describe the geometry of structure at a given cell is given by 

S = (A 2 ,A 3 , A 1 / ,A^,A 3 '). (7) 



3. Geometry with A -space distributions 

In Equations \5\ and [6] are crude relations that we expect the X' parameters to satisfy in the 
vicinity of specific large-scale structures. This can be made a bit more quantitative - filaments and 
walls will have simple 'fingerprints' in X'^-X^-X^ space (hereafter, A-space) that we can look for in 
real data. 
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3.1. Three-dimensional fingerprints 

In Fig. [H we plot the A-space projections for point distributions drawn from a two-dimensional 
(left), one-dimensional (centre), and three-dimensional gaussian function (right). These are models 
of a filament, wall, and clump, respectively, and are each placed in a 200 h~ l Mpc box. Smoothing 
will be done on a scale larger than the width of the structure, so the mock structures approximate 
a line, a plane, and a point. We confirmed that the fingerprints are coordinate-independent by 
rotating the mock structures to random orientations. 

Fig-EDplots the distribution of the A parameters evaluated at the position of each galaxy, each 
interpolated on a 128 x 128 x 128 grid. The points in this figure are colour-coded by local density 
(blue, green, yellow, orange, and red in order of increasing density). Note that the definition, 
A{ < Aj < A3, entirely excludes four octants in A -space, as well as portions of the remaining octants 
(delimited by the dashed lines in Fig.Q]). In this a nd all subsequent multidimensional in terpolations, 



we use a third-order polynomial routine given in lPress. Flannery. & Teukolskyl dl986). 



In all of the A -space distributions, there is a tight grouping of blue points near the origin made 
up of background objects well away from the mock structures; the Hessian matrix of a uniform 
density field is all zeros. The remaining points trace the 'fingerprint' of the wall, filament, or 
clump. Each fingerprint can be understood by considering the curvature properties of a gaussian 
function. A gaussian function has curvature given by 

where fix) is a one-dimensional gaussian function with width, o. At x = 0, the right-hand side 
of Equation [8] is equal to — p(x)/o 2 . Moving to larger values of \x\, the value of this function 
increases, reaching a maximum at |jc| = \/3<7, and then switching sign at \x\ = o. Beyond this, 
the curvature of a gaussian asymptotically approaches zero. The fingerprint of the wall follows 
naturally from this simple analysis. At all points near the wall centre (the densest regions, in red), 
the axis of structure will be either aligned with or orthogonal to the x coordinate axis - the density 
field has no y or z dependence so the curvature of the density field along those axes will always 
be zero. This means that will correspond to either X[ (the smallest eigenvalue) or A3 (the 
largest eigenvalue), depending upon the sign of the curvature along the x direction (which follows 
Equation [8] after smoothing). Therefore, for |x| < o 

A, = 0<O (9) 
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and for |x| > o, 

= = ^ = 

This explains the origin of the half-cross shape in the centre panel of Fig. [TJ 

For the clump, the direction of least curvature (and, therefore, the axis of structure) will be 
either parallel or perpendicular to the radial vector from the clump centre. Along this radial vector, 
the density field will follow the curvature properties of a gaussian function, corresponding to A3 
in the bottom-right panel of Fig. [Q Orthogonal to this vector, the density field will always have 
negative curvature, meaning that the smallest two eigenvalues (A{ and X^, see the top-right panel) 
will always be negative and, by symmetry, equal in magnitude. All three eigenvalues increase with 
distance from the clump centre (in red) until the curvature of the gaussian function maximizes at 
r = \/3(7. Beyond this, A3 drops and the fingerprint approaches the origin. 

Finally, in the mock filament (left column of Fig. Q]), there is no z-dependence in the density 
field , so one of the A -values will always be zero and the fingerprint will be restricted to one of the 
A -space coordinate planes. Near the filament centre (in red), the curvature is negative orthogonal to 
the filament axis. This means that the two non-zero eigenvalues are negative and the fingerprint lies 
in the A^-Aj plane. Both of these eigenvalues will increase with distance from the filament centre, 
with the one corresponding to the radial eigenvector following the curvature of the one-dimensional 
gaussian function (Equation [8]). At a cylindrical radius, R = y/3o, the curvature changes sign and 
this eigenvalue becomes positive (and, therefore, becomes A3), after which the fingerprint remains 
on the Ag-Aj' plane and behaves like a projection of the clump fingerprint. 



3.2. Wall and filament proportions 



Neither the real data nor the cosmological simulations will have such a simple structure as the 
idealized models described the last section. To improve the realism of the tests, we will give the 
walls and filaments a finite extent, as well as place multiple structures within a single box. 



The first set of tests (see Bond 2008) varied the length of the filament, ranging from six times 
to twice the smoothing length. As the filament size was reduced, the filament began to resemble 
a clump in A-space. Filaments also have a wide range in width, as the very largest structures in 
the universe (e.g., the Sloan Great Wall. iGott et al.ll2005|) have thicknesses of ~ 50 hr x Mpc in 
redshift space, much greater than the widths of 'typical' filaments visible in redshift surveys. We 
generated sample filaments to test this width dependence, varying it between half and twice the 
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smoothing length. We found that, even for very long filaments, the A-space distributions show few 
of the characteristic 'filament' features unless the structure is narrower than the smoothing length. 
Thus, single-scale smoothing is insensitive to structures much wider than the smoothing length, 
suggesting that it will be useful to measure the scale-dependence of filamentarity with a series of 
measurements on different smoothing scales. 

We have also considered the impact of nearby structures on the filament fingerprint. We find 
that fingerprint remains largely unchanged so long as the density of filaments is not so high that 
nearby structures are within one smoothing length of one another. In this limit, a smoothing length 
less than the average separation between structures is needed to distinguish them. 

We conducted similar tests for walls, varying their side lengths, widths, and number within a 
box. Points in the vicinity of a wall's edge will have A -parameters similar to those of a filament 
and parameters for points near the corners will be like those of the clump. As with the filaments, a 
great deal of scatter was introduced when the structures were made very wide or very numerous. 

Although the morphology of the A-space distributions can vary with a structure's width, 
length, or overdensity, there are discriminating features for each type of structure. Filaments can 
most easily be distinguished from clumps by a concentration around the A3 = axis in the A3— A{ 
projection, while walls can be distinguished from filaments by a concentration on the X' 2 = axis in 
the Aj-Aj' projection. In §[71 these facts will be used to probe structure in cosmological simulations 
and SDSS galaxy data. 



4. The simulation 

Before we present the A -space projections of the observed galaxy distribution, it is useful to 
explore the predictions of the standard cosmological model. A cosmological simulation allows us 
to develop our techniques on a data set not subject to redshift distortions, sparse sampling, and 
complicated window functions. Furthermore, we can simulate these effects in mock catalogues, 
providing a useful standard for comparison to the real SDSS data. 

We use a series of cosmological CDM /Y-bo dy simulations with Q, m = 0.29, = 0.71, Ofc = 



0.85, and h = Hq/(\QQ km s 1 Mpc *)= 0.69 dSpergel et al]l2007l) . The simulation is performed 



within a 200 h 1 Mpc box with 512 3 particles, each with mass, m p = 4.77 x 10 9 hT 1 M , and 
is run on a 92-node Beowulf cluster located at Princeton University. The positions and velocities 
of dark matter particles are output at six redshifts - z = 0, 0.3, 0.5, 1, 2, and 3. Because of 
memory restrictions, we split the simulation box into a 3 x 3 x 3 lattice in order to run a halo- 
finder and identify galaxies. Each subbox was padded with the particle data within 10 hr x Mpc of 
its boundaries. The padding size is an order of magnitude larger than the largest dark matter haloes, 
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so edge effects are negligible. 

We identified dark matter haloes with the HOP algorithm ( Eisenstein & Hutlll998 . hereafter EH98), 
which searches for peaks in the density field by moving from particle to particle until it reaches one that is 
in a denser environment than any of its neighbors. Densities are computed with an adaptive kernel with 
a length scale set by the Nd eas nearest particles, and the 'hopping' can occur to any of the A/hop nearest 
particles. Groups are then pruned, merged or disbanded depending on their peak and outer density thresholds 
(see EH98 for more details). Following EH98, we select A^ p = 16, A'dens = 64, N meige = 4, 5 pea k = 240, 
Caddie = 200, and 5 oute r = 80. In addition, we keep only haloes with N > 8 particles. These p arameters give 
halo distributions similar to those of a friends-of-friends algorithm Jlfuchra & Geller 1982b with a linking 
length of 0.2. We note that our results will be insensitive to the choice of the minimum halo mass as the 
faintest galaxies we will be studying (see § 16.11 ) have a mean halo mass of ~ 10 M Q , corresponding to 
N ~ 200 particles. 



5. Gaussian random fields 

In both the cosmological simulations and observed galaxy distribution, we wish to identify the mor- 
phological features that are unique to non-linear structure evolution; that is, features which are unlikely to 
appear in a gaussian random field. If we wish to reveal these features in A -space, then we should compare 
the Hessian parameters in observed or simulated data to those in a gaussian random field with the same 



power spectrum. For the non-line ar power spect r um, w e use the prescription of ISmith et all (12003!) . with the 



concordance model parameters of ISpergel et al.1 (12007). 



The left column of Fig. |2] shows a slice from the z = 99 initial conditions for the dark matter distribu- 
tion in a cosmological simulation, smoothed on a 5 hT x Mpc scale. In the top panel we plot the smoothed 
density field, while the smoothed X[ and X' 2 maps are given in the centre and bottom panel, respectively. The 
grayscale map of X[ gives a strong impression of filamentary structure. Note, however, that the X parameters 
are essentially edge-finders and both overdensities (whose edges are apparent in the X[ map) and underden- 
sities (X' 2 map) will have boundaries. In a Gaussian random field, the overdensities and underdensities are 
statistically identical, so the maps are morphologically similar. 

On the righthand side of Fig. |2 we show the same grayscale maps for the z = dark matter distribution 
in a cosmological simulation (see §0]), but with the X parameters (centre and bottom panels) normalized to 
local densitjQ Non-linear growth of structure amplifies overdense edges and the low-contrast structures 
tend to become washed out in the raw A -maps. After density normalization, however, the edges show a 
similar morphology to those in the X[ and X^ maps of the gaussian random field, suggesting that the outline 
of the network of filaments and walls is already in place in the initial conditions (as was first suggested by, 



2 Note that the gaussian random fields used in the left columns of Figs. [2] and [3] are not the seed fields for the 
simulations use in the right columns. 
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Bond et al. 1996). Also included in the bottom two rows of Fig. [2] are bars indicating the direction of the 
local axis of structure. In the simulations, the axis of structure is nearly perfectly aligned with the X[ edges 
in even the weakest strands, while in the gaussian random fields, there is no obvious correlation between the 
axis of structure and the orientation of the edges. The direction of the axis of structure does not appear to be 
correlated on much more than a smoothing length in the gaussian random field, while the simulations have 
structures that are aligned across the entire 200 h~ x Mpc box! 

6. Geometry of dark matter and mock galaxy distributions 

The two-point correlation function and power spectrum already provide a direct measure of the two- 
point properties of the large-scale distribution of matter, so we wish to extract from A -space only that 
information which cannot be obtained from two-point statistics. In the left column of Fig. [3l we show 
the A -space distribution (smoothed on a 10 h~ l Mpc scale) of a three-dimensional gaussian random field. 
The contours illustrate A -space distributions for points in the bottom 50 per cent of the real-space density 
distribution (black), between 50 and 75 per cent (blue), and the upper 25 per cent (red). The gaussian random 
field was sampled in a (200 Mpc) 3 three-dimensional grid in proportion to the local overdensity, which was 
normalized to match the amplitude of the ACDM non-linear power spectrum at each smoothing scale. No 
points were placed in grid cells having 8 < — 1 (< 0.1% of the volume in all cases presented here). Within a 
given panel of Fig.|3l there is little variation in the size or shape of the contours as a function of overdensity, 
indicating no density dependence of structural morphology. The plot is generated from only one realization 
of a gaussian random field, but finite volume effects are small at this smoothing scale. 

In the large-scale distribution of matter, structure grows due to the gravitational effects of the total 
matter distribution. As such, we compute A -space distributions at the positions of dark matter particles (not 
grid cells) in the z = 0, dark-matter-only simulation described in § |4j these are shown in the right column of 
Fig. [3] At high densities (red), the dark matter contours are extended because of the presence of overdense 
clumps. The low-density contours, however, are contracted relative to those of the gaussian random field. 
The largest underdensities in the gaussian random field have a positive curvature (and, therefore, positive A- 
values) comparable to the negative curvature of the overdensities, while dark matter underdensities are voids 
with a very small spread in density and curvature. If non-linear growth of structure leads to the formation 
of filaments and/or walls in the dark matter distribution, there should be an excess in the regions of A -space 
occupied by these fingerprints. 

In order to bring out these regions, Fig. 0] shows the A -space distributions of dark matter after sub- 
tracting a gaussian random field with the same power spectrum (i.e. the difference between the columns 
in Fig. [3]>. Here, grayscale density maps replace contours and there is no colour-coding for real-space den- 
sity, but only points with density greater than the mean are included. A comparison of these plots with 
the fingerprints in Fig. Q] will help us determine which structures dominate the dark matter distribution at 
each smoothing scale. Note that the axes of each plot have been scaled in proportion to the amplitude of 
the dimensionless power spectrum at the appropriate smoothing scale, allowing for a direct comparison of 
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structural morphology between scales. 

In order to make this comparison more quantitative, we will define a set of structural parameters, 

4 = i-^> ("> 

where the mean is computed over the region of A -space having A < and Aj < 0. For example, R\ 2 quantifies 
the mean axis ratio for structures that are negatively curved along at least two principal axes of the Hessian 
ellipse, allowing one to distinguish walls from filaments. A small value of R\ 2 indicates that the axes have 
similar curvature and are more like filaments, while R\ 2 ~ 1 indicates the dominance of wall-like structures. 
Similarly, can be used to distinguish filaments from clumps in regions with negative curvature along all 
three principal axes. 

As shown in the righthand column of Fig. [TJ the densest parts of clumps lie on the lower boundary of 
A-space (dashed line) in all three projections. Filaments concentrate on the lower boundary of the A^-A{ 
projection, as do clumps, but filaments also concentrate along the x axis of the other two projections. In 
Fig. 01 we see that at both 15 and 10 h~ l Mpc scales, the points in Fig. @] appear to be shifted up toward 
the x axis relative to the clump fingerprint. The distribution of points in the A3 vs. A[ projection is broad 
and extends to very negative values of X[, consistent with the presence of filaments in addition to clumps. 
At z = 0, it appears that filaments are the dominant over clumps at 15 hT 1 Mpc scales (R 23 = 0.605) and 
10 h~ x Mpc scales (R^ = 0.602), with clumps becoming more prominent at 5 h~ x Mpc smoothing scales 
(/?23 = 0.554). This evolution is most apparent in the highest contrast structures (small A2). 

There is little evidence for walls at z = at the smoothing scales shown. Walls manifest themselves as 
an excess of points near both X 2 = and A3 = 0. Although they are difficult to distinguish from filaments 
in the A3 vs. A{ projection, there would be a clear overdensity of points along the x axis of the top plots of 
Fig. [3] and near the origin of the bottom plots if walls were dominant features. The mean eigenvalue ratios 
suggest that structures become less wall-like at 5 h~ l Mpc smoothing scales, with R\ 2 = 0.370, as compared 
to R\ 2 = 0.485 at 15 hr l Mpc smoothing scales. 

In Figs. [5]and[6l we present the A-space distribution of structures at z = 1 and z = 3, respectively, 
smoothing on the same three comoving length scales as in Fig. |4](z = 0). The morphology of overdense 
structures in the dark matter distribution changes significantly at these redshifts. In all three A-space pro- 
jections, the distribution is concentrated closer to the x axis at higher redshift, meaning that the larger 
eigenvalues decrease more rapidly than the smaller ones and overdensities are approaching spherical sym- 
metry. 

At I = 15 h~ l Mpc, all of the plots are more sparsely populated at high redshift, indicating that the 
A-space distribution is closer to that of a gaussian random field. This is especially true of A3 vs. A^, which 
also shows no dramatic change in shape out to z = 1. In the A3 vs. A{, / = 15 h~ l Mpc projection, the most 
densely populated regions retain their shape and relative contrast out to the highest redshift (/?| 3 = 0.794 at 
z = 3 as compared to 7?j 3 = 0.764 at z = 0), but some of the sparse features in the bottom left corner fade 
by z ~ 1- Finally, the upper left panel, X 2 vs. A{, undergoes strong evolution, with the tails at very negative 
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eigenvalues disappearing and the bulk of the distribution moving to the A^ = axis. This evolution is also 
demonstrated in R\ 2 , which drops from 0.599 at z = 3 to 0.485 at z = 0. 

In Fig. [T] we see that walls are most easily distinguishable from other structures in A^ vs. A[, with 
points concentrated along X 2 = 0. Taken together, these facts are suggestive of wall-to-filament evolution 
at a / = 15 hT x Mpc comoving smoothing scale. The rapid fading of the A3 vs. X 2 map is also explained, 
as walls are concentrated toward the origin in these two eigenvalues and would be indistinguishable from a 
gaussian random field. The signature of walls at high redshift are apparent in these plots, but they are of low 
contrast compared to the structures seen at smaller scales and lower redshift. 

The two smaller smoothing lengths, I = 10 h~ x Mpc and / = 5 hT x Mpc, exhibit a filament-to-clump 
evolution. In A3 vs. A£, the distributions at both length scales exhibit a tail along the dashed boundary that 
grows toward z = 0. This tail is populated by the points near the dense centres of clumps. A similar tail 
exists in A3 vs. A(, but here the growth of individual filamentary features is also evident, with a concentration 
near the A3 = axis. The eigenvalue ratios, R 23 = 0.647 at z = 3 and R 23 = 0.554 at z = 0, suggest filament 
dominance at / = 5 hT y Mpc across the entire redshift range, but with an evolution towards more "clump- 
like" filaments. 

It is interesting to compare the z = 1 projections for / = 10 h~ l Mpc to the z = 0, 1 = 15 h~ l Mpc pro- 
jections. Their similarity suggests that the two length scales are in similar stages of non-linear evolution, 
but at different times. This conclusion is supported by the eigenvalue ratios, with / = 15 h~ l Mpc having 
R^ 2 = 0.485 at z = and I = 10 hr x Mpc having R\ 2 = 0.497 at z = 1. That the smaller length scales develop 
filamentary structure earlier is an indication of bottom-up structure formation; however, at a given comoving 
scale, struc t ure ev olves from walls to filaments to clumps, consistent with the ellipsoidal collapse picture of 
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6.1. Mock catalogues in A -space 



In the SDSS data, we will be working with the observed galaxy distribution, so it is useful to generate 
A-space projections for mock galaxies, generated by the methods described in § [4j There are many fewer 
galaxies than dark matter particles in the box volume, so the signal will be much weaker and shot noise 
may be non-negligible at / = 5-10 \T X Mpc scales. More importantly, there may be interesting differences 
between the large-scale structures appearing in the dark matter and galaxy distributions. In order to pop- 
ulate the dark matter haloe s with galaxies, we must know the relevant halo occupation distribution (HOD, 
Berlind & Weinbergll2002r) : the probability, P(N\M), that a halo of mass M will contain N galaxies. 



Here we use the parametrization of IZheng. Coil. & Zehavil ((2007), in which the probability of a halo 
of mass M containing a central galaxy is given by: 



(AT. 



1 + erf 



logM - logM n 



(12) 



Olog M 

Each central galaxy may be surrounded by one or more satellite galaxies, the mean number of which is given 
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by, 



(AT,) = i 1+erf 



( 



logM - logM, 

Clog M 



) 



( 



M-M 
Mi 



) 



a 



(13) 



The number of satellite galaxies in a given halo is drawn from a Poisson distribution and has no upper limit. 
The free parameters, M m j n , O\ og m, M)> Mi, and a have been derived from fits to SDSS galaxy catalogues. 
Table 1 of Zheng et al. 2007 provides the best-fitting HOD parameters for a series of absolute r-band 
magnitude cutoffs. 

Once drawn from the relevant distributions, central galaxies are given the mean position and peculiar 
velocity of their associated dark matter halo, while satellite galaxies are each given the position and velocity 
of a randomly selected dark matter particle within the halo. When working with the SDSS data, we will con- 
struct a series of volume-limited samples, so the mock catalogues must be prepared to mimic the associated 
selection function. 

One important difference between the redshift-space galaxy distribution and the dark matter distri- 
bution is the presence of fingers-of-god. These features, created by the peculiar velocities of galaxies in 
clusters, appear as very narrow filaments along the line of sight with widths ~ 1 hr x Mpc and lengths up 
to ~ 10 /r 1 Mpc. The mock catalogues must contain a model of these redshift distortions if they are to 
properly mimic the SDSS data. The 'distorted' positions of galaxies are computed, 



where r is the position in real space, v is the galaxy's peculiar velocity, and r is the unit radial vector between 
the observer and the galaxy. 

In Fig.|7J we plot the A-space projections of an M r < —20.5 mock galaxy sample projected into redshift 
space. The redshift distortions are applied from a vantage point with coordinates {-300,0,100} h^ 1 Mpc rel- 
ative to the origin of the box, selected to closely resemble the data samples constructed in the next section. 
The morphological characteristics of the A-space distributions for / = 15 and 10 h~ x Mpc are very similar 
to those of Fig. |U though several features are different (for example, the smearing of the previously noted 
clump signature in the bottom left corner of A^ vs. X[). The changes are most noticeable at / = 5 h~ l Mpc, 
where all of the increase by ~ 0.05 in the mock catalogs and the high-contrast features all but disap- 
pear. These high-contrast features are primarily due to clusters, which become lower-contrast fingers-of-god 
when projected into redshift space. 

With the construction of the mock catalogues and the development of the theoretical framework pre- 
sented in § |2j we now have all of the tools necessary to compare the real galaxy distribution to that predicted 
by the standard cosmological model. In the next section, we will present our galaxy sample, drawn from the 
SDSS, and use the A -space distribution to make inferences about real large-scale structure. 



r-v ^ 



(14) 
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7. Filaments and walls in the SDSS galaxy distribution 



The Sloan Digital Sky Surve y has imaged a c 



to 10000 A, to a depth of r ~ 22.5 (lAbazajian et al. 



uarter of the sky in five wavebands, ranging from 3000 



2009). As of Data Release 7, spectra have been taken of 



930,000 galaxies, covering 8032 deg 2 and extending to Petrosian r ~ 17.7 (lStrauss et al.ll2002r) . Galaxy 



redshifts are typically accurate to ~ 30 km s~ 1 , making the survey ideal for studies of large-scale struc- 
ture. For this study, we need a portion of sky with relatively few coverage gaps to minimize the effect of 
the window function on the A-space distributions. With this in mind, we construct three volume-limited 
subsamples from the north ern portion (8 < a < 16 h and 25 < 8 < 60) of the NYU Value-Added Galaxy 
Catalogue (NYU-VAGC, blanton et allboftd , through DR6), the first 80 x 80 x 220 (/T 1 Mpc) 3 in size 
with M r < -20 (Mr20), another 140 x 140 x 340 (/r 1 Mpc) 3 in size with M r < -20.5 (Mr205), and finally 
170 x 170x400 (/T 1 Mpc) 3 in size with M r < —21 (Mr21). Absolute magnitudes were computed with 
kcorrect using SDSS Petrosian magnitudes shifted to z = 0.1 (and using h = 1). The sample properties of 
the galaxy distribution are given in TableQ]and the galaxy distributions are plotted in Fig. [8] 

All boxes are padded with 30 h~ l Mpc of empty space and binned on a 128 x 128 x 256 grid, resulting 
in grid spacings of 1.09, 1.56, and 1.80 hT l Mpc in the Mr20, Mr205, and Mr21 subsamples, respectively. 
To correct for edge effects, we define a window function on the same grid such that, 



Wi 



1 if ieS 
if i^S 



(15) 



where 5 is the set of grid points within the SDSS survey volume and i is a grid cell within the subsample 
volume. The above window function was computed by smoothing the sample boxes on a 20 h~ 1 Mpc scale 
and including in the window only those regions with densities more than 4 standard deviations below the 
mean density in the box. This is a crude approximation, but sample boxes were selected to lie in regions 
with few SDSS coverage gaps, so only a small fraction of grid cells are excluded (primarily near the box 
boundaries). 

Once the window function and density fields have been computed on the grid, we smooth both with 
gaussian kernels on the scales given in Table [Q giving W(x) and p(x), respectively. The corrected density 
field is then 



p\x) 



W(x)' 



(16) 



The first and second derivatives of the density field have zero mean on large scales (if the cosmological 
principle holds) and the kernels used to compute them will integrate to zero. As such, there is no need to 
renormalize the smoothed derivative fields near the padded regions. Rather, we need only subtract the signal 
from the edge artificially introduced by the padding. That is, 

H' li =H ii -pB lh (17) 
where is the matrix of smoothed second partial derivatives of the window function. 
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For each of the data samples listed in Table [TJ we constructed a mock galaxy catalogue that matched 
the box dimensions, number density, and padding of the SDSS data, incorporating redshift distortions ac- 
cordings to the prescription described in § 16. 1[ The processing on the mock galaxy catalogues and the real 
data were identical and they should be directly comparable, with the caveat that the mock catalogues were 
taken from a 200 x 200 x 200( h~ x Mpc) 3 box periodic box, and therefore had to be repeated along the z 
axis in order to match the dimensions of the data subsamples. As such, finite-volume uncertainties will be 
larger in the simulations. 

In Fig. |9j we compare the two-dimensional X[ maps between 10 h~ x Mpc thick slices from the M r < 
—20.5 subsample, its corresponding mock catalogue, and the dark matter distribution. These data represent 
~ 3 per cent of the M r < —20.5 subsample and the normalization and range of the grayscales are identical 
between the two columns of each figure. The mean properties of each slice are superficially similar, although 
the slice from the mock catalogue has a large filament running through its centre, dominating the visual 
impression. In all maps, the axis of structure is oriented along edges in \[, unlike the gaussian random field 
shown in Fig. [2] 

The A-space projections of the SDSS data are displayed in Fig. [l0]for / = 20 h~ x Mpc (Mr21), / = 
10 hT x Mpc (Mr205), and / = 5 hT x Mpc (Mr21). For comparison, we present the corresponding A-space 
distributions for the redshift space mock catalogues in Figs.Qj] At small scales, the value of A{/A^ at the 
positions of galaxies is generally near unity (diagonal lines) and the distribution morphologies are similar to 
those in the mock catalogues. Shot noise is more severe at these scales than at larger scales, so the similarity 
could be in part due to the extra power added by Poisson fluctuations. On 10 lr x Mpc scales, there are some 
indications that the real galaxy distribution may be more evolved than the simulated one, with all R x being 
smaller in the data by ~ 0.05. The A-space distributions for the data appear to have more high-contrast 
structures, including a tail of clump-like objects with near-unity values of A^/A^. If real, the discrepancy 
could suggest that the simulation is underestimating o&. 

Finite- volume effects are apparent at the larger smoothing scales for all of the samples, where the Re- 
values are sometimes larger and sometimes smaller than their counterparts in the mock catalogues. At the 
largest smoothing scales, there is clear evidence for non-linearity and the A^-A{ projection shows evidence 
for both filaments and walls. Unlike the distributions at / = 5 hT x Mpc and / = 10 h Mpc, the data do not 
appear to be much more clump-like than the mock catalogues (R^ =0.61 in the data and R^ 3 = 0.54 in the 
mock catalogues). However, there is some suggestion that the data are less wall-like on the largest scales 
(R\ 2 is smaller in the data by ~ 0.08). This would also be consistent with a more evolved population, though 
finite-volume effects and uncertainties in the halo occupation distribution make it difficult to say for sure. 
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8. Discussion 
8.1. A matter of scale 

Mention of large-scale structure often brings to mind images of narrow, interconnected filaments with 
clusters at the interstices weaving their way between vast cosmic voids. Bond et al. 1996 dubbed this mental 
image the 'cosmic web', but we have avoided using that phrase in this paper. The components of a spider 
web are strands, which come together at nodes to make a variety of complex structures on larger scales. 
The same could be said for filaments in the dark matter distribution, but that is not a complete picture. We 
also know that these filaments are made up of bound dark matter haloes. These small-scale filaments can 
themselves be the building blocks for large-scale ones. 

Some recent attempts to quantify and identify individual filamentary structures have n ot taken the mul- 



tiscale nature of struc ture into account. For example, minimal spanning tree algorithms (e.g. jBarrow. Bhavsar. & Sonoda 



1985: IColbergl 120071 ^ attempt to select filaments in the same manner as clusters, assuming that a step down 
up in linking length (that is, down in density contrast) will allow identification o f the strands and full recon- 
struct ion of the web. Adaptive smoothing kernels and tesselation methods (e.g., Schaap & van de Weygaertl 



booch enhance high-density regions in much the same way that simulations of galaxy formation adapt the 
mesh to treat star-forming regions, resulting in a nebulous combinati on of structure on small sc ales and large 



scales. There have also been attempts to use image segmentation (lAragon-Calvo et alj|2007f) to select the 
most 'prominent' large-scale structures across a multitude of scales, and then remove them before looking 
for the next ones. Although this procedure makes sense on a given scale, removal of structures on one scale 
can cause one to miss structures on another. 

We demonstrate the problem with this last approach using a 15 h~ [ Mpc deep slice from the SDSS 
Mr20 subsample (Fig. [12]). At smoothing scales of both 10 h~ [ Mpc (upper right panel) and 3 h~ l Mpc (lower 
right panel), the axis of structure appears to be aligned with the A[ ridges, suggesting that matter is condensed 
into filaments/walls on both scales (see Fig. [2]). Furthermore, the filaments/walls found on the smaller scale 
appear to be components of those on larger scales. If the structures found on the 10 hT x Mpc scale were 
removed before any of the smaller structures were identified, information would be lost. Conversely, if the 
smaller structures were all identified and removed first, there would be little left to find the larger ones. 

In this paper, we have presented results on a number of smoothing scales and demonstrated how struc- 
ture changed across those scales. Filaments appear to be present on scales ranging from at least 5 h~ x Mpc to 
25 h~ x Mpc and are the dominant structures throughout most of that range. In Paper II, we extend these 
ideas by presenting a method to identify and characterize individual filaments. 



8.2. Walls, walls, everywhere 



There is evidence for a filament-to-clump progression as one smooths on progressively smaller scales 
in the z = A-space projections in both simulations and real data. Since smaller-scale overdensities collapse 
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before larger ones, they should be in later stages of growth. The collapse of a large-scale overdensity might 
be expected to follow such a progression even if it were composed of smaller-scale condensed objects (such 
as dark matter haloes). The coherence of these large-scale filaments and walls would presumably depend 
upon the degree of fragmentation on intermediate scales. 

One means of testing this picture is to look for walls in the very early stages of non-linear evolution. 
The z = A -space projections of the dark matter distribution show no convincing evidence for wall-like 
structures at scales as large as 15 h~ x Mpc, but they may still exist on larger scales or earlier times. The 
former is difficult to explore in our 200 h~~ 1 Mpc simulation box, but we do have simulation outputs out to 
z = 3. In Fig. [13] X[ maps for a 15 h~ l Mpc smoothing length are plotted for a range of redshifts, along 
with bars indicating the axis of structure. Even at z = 3, it is clear that, unlike in a gaussian random field, 
the axis of structure aligns with the X[ edges (though it is impossible to distinguish filaments from walls in 
a two-dimensional projection). As discussed in §|6]and § 18.11 there appears to be a progression in both time 
and smoothing scale of wall-to-filament-to-clump. If this progression holds at z = 0, we might expect to see 
wall-like structures dominating on ~ 40-50 hr x Mpc scales. 

The panels in the left-hand column of Fig. [14] show a subsample of dark matter particles in the 
200 hT x Mpc simulation box satisfying X[ < / C ut<7i5(z)> where dj 2 5 (z) is the variance of the density field 
at a given redshift when smoothed on a 15 h~ l Mpc scale. Here, X[ is measured with 15 /i 1 Mpc smooth- 
ing and / cut is a dimensionless constant. The cut isolates very similar regions of space at each redshift, but 
the distribution of matter within those structures evolves a great deal with time. The regions isolated by the 
cut are very sheet-like, and at z = 3 the matter is smoothly distributed within the structures. As the evolution 
progresses, they collapse into filamentary and clump-like structures. 

Similarly, the right-hand column Fig. [14] shows particles in regions of small X[, but in a smaller box 
and on a smoothing scale of 5 h~~ 1 Mpc. Here, structures have begun to collapse into filaments at z = 3 and 
by z = are primarily concentrated in clumps. Despite the evolution of the structures within the regions 
selected by a cut in X[, the sizes and morphologies of the regions themselves change very little - they are 
still pancake-like at z = 0. This suggests that, even at low redshift, there should be evidence for wall-like 
structures in the galaxy and dark matter distributions on large scales, even if these walls are themselves 
made up of prominent filaments and clumps. As such, at all redshifts probed here, it can be the case that the 
majority of the mass in the universe is in walls on large scales, which are made up of filaments on smaller 
scales, which are made up of clumps on still smaller scales. 

8.3. Conclusions 

With the eigenvalues and eigenvectors of the Hessian matrix, we can determine the type and orientation 
of structures in a continuous two- or three-dimensional density field. The 'fingerprints' of clumps, walls, 
and filaments in a three-dimensional space of eigenvalues (A -space) provide templates for comparison to 
the large-scale matter distribution. In addition, as a result of non-linear growth of structre, the third Hessian 
eigenvector (which is used to define the 'axis of structure') aligns with the axis of filaments in both two and 
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three dimensions, retaining this alignment within approximately one smoothing length from the central axis. 

We drew three volume-limited subsamples from the northern portion of the SDSS spectroscopic survey 
(using the NYU-VAGC catalogue) and computed their Hessian parameter distributions on a range of length 
scales. These distributions were then directly compared to those found in a series of redshift-space mock 
galaxy catalogues generated from a cosmological simulation using the concordance cosmology. Results 
from this analysis include, 

• The SDSS galaxy distribution is different from a gaussian random field on all smoothing scales used 
here (up to 25 h^ 1 Mpc), showing evidence for some combination of clusters, filaments, and walls on 
all of these scales. 

• Filaments are the dominant structures on smoothing scales of ~ 10 h~ x Mpc to ~ 25 h~ l Mpc, though 
there is some evidence for a transition to wall-dominance at ~ 25 h~ x Mpc. At 5 h~ l Mpc, clumps 
begin to dominate for the highest-contrast structures, but filaments are still apparent in the A -space 
distributions. 

• In very negatively curved regions of the galaxy density field (as would be found near walls, filaments, 
and clusters), the direction of least curvature, indicated by the Hessian eigenvectors, points to other 
very negatively curved regions. In other words, the axis of structure aligns with minima in X[. This 
occurs in projected slices at all smoothing scales studied here, providing a signature of non-linear 
growth even when non-linearities are not apparent to the eye in the density field. 

• The A-space projections of the SDSS galaxy distribution are very similar to those expected in a 
ACDM universe with gaussian random phase initial conditions. There is some evidence that the data 
are clumpier than the simulations. 

• In projected slices, the A[ maps of the galaxy distribution are morphologically very similar to those 
of a gaussian random field. Since filaments and walls are oriented along minima in %[, this suggests 
that the outline for the filament network can be determined from the initial conditions. 

To complement the results presented above for the SDSS galaxy distribution, we generated A pa- 
rameters at six redshifts in a ACDM cosmological N-body simulation, tracing the evolution of the matter 
distribution from z = 3 to z = 0. On comoving smoothing scales of 5 h~ x Mpc, filament-dominance at 
z = 3 gives way to clump-dominance at z = 0. At z = 3 and on comoving scales of 15 h~ l Mpc, the dark 
matter distribution is dominated by low-contrast walls, despite being only barely distinguishable from a 
gaussian random field in A -space. Structure on this smoothing scale becomes filament-dominated by z = 0. 
Finally, in projected slices, the axis of structure follows the minima in X[ on smoothing scales as large as 
15 /z 1 Mpc and epochs as early as z = 3. The same alignment is not seen in gaussian random fields. 



The A -space distributions of the galaxy de nsity field are only a first ste p towards a more complete 
description of large-scale structure. In Paper II (|Bond. Strauss. & Cenl l201fj) . we describe a method for 
finding individual filamentary structures on a given scale, and then quantitatively evaluate the connection 
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between filaments and galaxy clusters. Eventually, we hope to both extend our morphological studies to 
larger redshifts and to characterize the galaxy populations in structures on a range of scales. We will compare 
the SDSS A-space distributions to those found in hydrodynamic simulations, test the dependence of the X- 
space distributions on t he cosmological model and galaxy HOD, look for a dependence of galaxy properties 
on th e local geome try (|Park et al.ll2007h . and compare the galaxy A-space distribution to that of the hot gas 
re.g. lphillipsll2003h . 
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Tablel. SDSS Subsamples 
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Fig. 1. — The structural 'fingerprints' for, from right-to-left, an isolated filament, a wall, and a 
cluster in an otherwise uniform background, displayed in three projections of A -space. The points 
are colour-coded by local density (blue, green, yellow, orange, and red in order of increasing 
density) and inaccessible regions are shaded light gray, with eigenvalues defined such that X[ < 
< A3. All structures are smoothed on a scale much larger than their width. 
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Fig. 2. — A slice (10 h~ Mpc deep) from the dark matter initial conditions of a cosmological 
simulation at z = 99 (left column) compared with the same slice at z = (right column). The top 
row shows the density fields smoothed on a 5 Mpc scale, while the middle and bottom rows 
show the corresponding X[ and X ! 2 maps, respectively. The A maps have been normalized to the 
local density for the z = dark matter distribution in order to emphasize the large-scale distribution 
of structures and provide an easy comparison to the initial conditions in the left column. Bars have 
been plotted over the A -maps to indicate the direction of the axis of structure at random points on 
the plane. The bar's length is proportional to the magnitude of X[ at that point. It is clear to the eye 
that the axis of structure aligns with the local filamentary structure in the dark matter distributions. 
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Fig. 3. — The A-space projections for a three-dimensional gaussian random field (left) and the 
simulated z = dark matter distribution (centre). Contours represent the distribution of A at the 
positions of dark matter particles (not grid cells) from a sampling of the fields. In the left two 
columns, contours are for different density subsamples, where black contours indicate the dis- 
tribution of points in underdense regions (first or second quartile), blue contours in regions of 
intermediate density (third quartile), and red contours in the most dense regions (upper quartile). 






Fig. 4. — The A-space projections for the z = dark matter distribution, after subtraction by the 
A -space distribution of a gaussian random field with the same power spectrum. Grayscale cells 
indicate those regions in which the dark matter distribution shows an excess over the gaussian 
random field and for which the local density is greater than the mean, with darker cells indicating 
more of an excess. Regions showing no excess or regions showing a deficit are left blank. The axes 
of each plot have been scaled in proportion to the amplitude of the dimensionless power spectrum 
at the appropriate smoothing scale. 








Fig. 5. — The same as Fig.0 but for the z = 1 dark matter distribution. 




Fig. 6. — The same as Fig. @] but for the z = 3 dark matter distribution. At 15 h Mpc and 
10 hT x Mpc scales, distributions appear more wall-like than those at z = 0, suggested by a 
shift towards the A2 = axis in the A2-A1 projection. Meanwhile, the A -space distribution at 
5 hT x Mpc scales is more filamentary, as suggested by a similar shift toward the A3 = axis in the 
bottom two projections. These latter projections are very similar to those for 15 /z 1 Mpc scales at 
z = 0. 




Fig. 7. — The A-space projections for a z = mock galaxy catalogue with M r < —20.5, including 
redshift distortions. 
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Fig. 8. — The three volume-limited samples taken from the SDSS VAGC large-scale structure 
sample, with galaxies placed at their comoving positions based on the concordance cosmology. 
The location of the Milky Way is r = (250, -20, 1 10) h~ x Mpc, r = (310, -20, 170) h~ l Mpc, and 
r = (400, -25,200) hr l Mpc in the Mr20, Mr205, and Mr21 samples, respectively, and the z axes 
are parallel to the Galactic North pole. 
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Fig. 9. — The two-dimensional X[ distributions in slices from the z = output of a ACDM 
cosmological simulation (dark matter and mock galaxy distributions, left and centre columns), 
and the SDSS galaxy distribution (right column), smoothed on a 15 hr l Mpc scale (top row), a 
10 hT x Mpc scale (middle row) and a 5 hr x Mpc scale (bottom row). Bars have been plotted 
over the maps to indicate the direction of the axis of structure at each point in space. In all cases, 
the bar's length is proportional to the magnitude of X[ at that point. The grayscales are identical 
between the three columns. In all three samples and on all three smoothing scales, the axis of 
structure follows the minima in X[. 






Fig. 10. — The A-space projections for the three volume-limited subsamples of galaxies in the 
VAGC galaxy catalogue, smoothed on 20, 10, and 5 h~ l Mpc scales. Grayscale maps are the same 
as in Fig.Sl 
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Fig. 11. — The same as Fig. [TOl but for the mock galaxy samples. Grayscales and axis limits are 
the same between the two figures. 
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Fig. 12. — A two-dimensional demonstration of the structural hierarchy in a 15 h~ Mpc deep 
slice from the SDSS Mr20 subsample. Galaxies are plotted in the upper left panel and X[ grayscale 
maps are shown in the other two panels, smoothed on a 10 hr x Mpc scale (upper right) and a 
3 h Mpc scale (lower right). Bars over the grayscale maps indicate the local direction of the 
axis of structure and their length is proportional to X[. Alignment of the axis of structure with 
the X[ minima suggests that the structures are real on both scales (see Fig. |2|). The filament- or 
wall-like structures seen on a 10 hr x Mpc scale are divided into smaller filaments or walls on the 
3 hr 1 Mpc scale. 






Fig. 13. — The X[ grayscale maps for a series of 10 h~ l Mpc deep slices from the dark matter 
distribution, very similar to those in Fig. [2j The smoothing scale is 15 hr x Mpc and maps are 
plotted for z = 0, 0.5, 1, and 3. Bars indicate the direction of the axis of structure at random points 
on the grid, and the bar length is proportional to the local value of X[. In all maps shown here, the 
grayscale is scaled to the minimum and maximum of X[ in the z = map. The axis of structure 
aligns with minima in X[ out to z = 3, unlike in gaussian random fields. 




Fig. 14. — Selected dark matter particles in a 200 h Mpc simulation box at z = and z = 3. 
In the left column, particles are selected by the local value of X[/o\s{z) after 15 hr x Mpc gaus- 
sian smoothing. This cut follows individual structures as they evolve. Most of the structures 
are sheets/walls at z = 3, but have collapsed into filaments by z = 0. In the right column, 
5 hr x Mpc smoothing is used in a 67 hr x Mpc box. On that scale, most of the structures are 
filamentary at z = 3, but have collapsed into clumps by z = 0. 



